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ABSTRACT 

A scenario in which cosmic rays (CRs) above 10 20 eV are produced by cosmological 
gamma-ray bursts (GRBs) is consistent with observations provided that deflections 
by the inter-galactic magnetic field (IGMF) delay and spread the arrival time of the 
CRs over > 50yr. The energy lost by the CRs as they propagate and interact with the 
microwave background is transformed by cascading into secondary GeV- TeV photons. 
We show that a significant fraction of these photons can arrive with delays much 
smaller than the CR delay if much of inter-galactic space is occupied by large-scale 
magnetic "voids", regions of size ^ 5Mpc and field weaker than 10 _15 G. Such voids 
might be expected, for example, in models where a weak primordial field is amplified 
in shocked, turbulent regions of the intergalactic medium during the formation of 
large-scale structure. For a field strength ~4x 10~ n G in the high field regions, the 
value required to account for observed galactic fields if the IGMF were frozen in the 
protogalactic plasma, the delay of CRs produced by a burst at a distance of lOOMpc 
is ~ lOOyr, and the fluence of secondary photons above lOGeV on hour-day time 
scales is /(> E) ~ 10 _6 £ , ^ e 1 v cm -2 . This fluence is close to the detection threshold of 
current high-energy 7-ray experiments. Detection of the delayed flux would support 
the GRB-CR association and would also provide information on the IGMF structure. 



Subject headings: cosmic rays — gamma rays: bursts — magnetic fields 



Introduction 



Recent gamma ray and cosmic ray observations give increasing evidence that the sources of 
gamma ray bursts (GRBs) and of cosmic rays (CRs) with energy E > 10 19 eV are cosmological 



(see Fishman fc Meegan 1995] for GRB observations review; Bird et al. 1994, Yoshida et al. 



|1995 , Waxman 1995b for CRs). The sources of both phenomena, however, remain unknown. In 
particular, most of the CR sources discussed so far have difficulties in accelerating CRs up to the 



highest observed energies (e.g., Cronin 1992). Although the source of GRBs is unknown, their 
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observational characteristics impose strong constraints on the physical conditions in the 7-ray 
emitting region (Piran 1994, Meszaros 1995| ), which imply that protons may be accelerated in this 
region to energies 10 20 — 10 21 eV (Waxman 1995a, Vietri 19"95| ) . In addition, the average rate (over 
volume and time) at which energy is emitted as 7-rays by GRBs and in CRs above 10 19 eV in the 
cosmological scenario is, remarkably, comparable ( Waxman 1995a| ,b). These two facts suggest that 
GRBs and high-energy CRs may have a common origin. 

An essential ingredient of the GRB model for CRs is the time delay due to intergalactic 
magnetic fields. The energy of the most energetic CR detected by the Fly's Eye experiment is in 
excess of 2 x 10 20 eV ( Bird et al. 1994|) , and that of the most energetic AGASA event is above 
10 20 eV ( Yoshida et al. 1995] ). On a cosmological scale, the distance traveled by such energetic 



particles is small: < lOOMpc (50Mpc) for the AGASA (Fly's Eye) event (e.g., |Cronin 1992| ). Thus, 
the detection of these events over a ~ 5yr period can be reconciled with the rate of nearby GRBs, 
~ 1 per 50 yr in the field of view of the CR experiments out to lOOMpc in a standard cosmological 
scenario (e.g., |Cohen fc Piran 1994 ), only if there is a large dispersion, > 50yr, in the arrival time 
of protons produced in a single burst (Note, that this implies that if a direct correlation between 
high energy CR events and GRBs, as recently suggested by Milgrom & Usov (1995), is observed 
on a ~ lOyr time scale, it would be strong evidence against a cosmological GRB hypothesis). The 
required dispersion may result from deflections of CR protons by the inter-galactic magnetic field 



([Waxman 1995a|) . 



The inter-galactic magnetic field (IGMF) has not been detected so far. Faraday-rotation 
measures set an upper limit of ~ 10 -9 G for a field with IMpc correlation length (see Kronberg 1994 
for review). Other methods have recently been proposed to probe fields in the range 

10-io_i -20 G 



(e.g., Plaga 1995 , Lee et al. 1995] , |Kosowsky Sz Loeb 1996 ). Theoretical considerations regarding 
the existence and strength of the IGMF are related to the formation of the observed ^G fields in 
galaxies. Recent studies suggest that a galactic dynamo cannot produce the observed large-scale 
fields in galactic disks ( |Kulsrud fc Anderson 1992j ) and that one must turn to alternative 
mechanisms, which typically rely on a pre-existing field. Galactic fields might be created, for 
example, by compression of much weaker fields in collapsing protogalactic regions. This mechanism 



requires a protogalactic field of strength 10 -10 G and correlation length of order IMpc. 
Such fields could be primordial, in which case they would likely permeate all intergalactic space. 
However, this need not be the case. For example, such fields could be generated from a much 
weaker primordial field, ~ 10~ 20 G, due to the turbulence induced in the formation of large scale 
structure in the universe ( [Kulsrud et al. 1995| ). In this picture, the IGMF would "trace the 
mass", with high 10~ n -10 -10 G fields in the high density (proto-galactic) regions of the large scale 
structure and very low fields in the intervening voids. 

Most of the energy lost by the CRs as they propagate and interact with the microwave 



background is transformed by cascading into secondary GeV-TeV photons (e.g., Aharonian et al.\ 
1994, Protheroe fc Johnson 1996| ). In this Letter we show that even though the CR time delay 
must be ^ 50yr, a significant fraction of the GeV-TeV cascade radiation can arrive with much 
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shorter delays, on the order of hours to days, provided that a large fraction of the inter-galactic 
space is occupied by magnetic "voids", regions of very low magnetic field (< 10~ 15 G). In §2 we 
present a qualitative discussion of the development of electro-magnetic cascades in the presence 
of an IGMF. The expected high energy photon flux is calculated using detailed Monte-Carlo 
simulations in §3. Implications for current and future high energy gamma-ray experiments are 
discussed in §4. 

2. Electro-magnetic Cascades and the Inter-galactic Magnetic Field 

As they propagate, nucleons with energy in excess of 10 20 eV lose energy due to the 
photo-production of pions in interactions with microwave background photons. Pion decay 
converts ~ 40% of the energy lost by nucleons to neutrinos, and the rest to photons, electrons and 
positrons. The energetic, ~ 10 19 eV, secondary photons, electrons and positrons further interact 
with the microwave background to form electromagnetic cascades: high energy photons interact 
with background photons and produce electron-positron pairs, which, in turn, lose their energy by 
inverse-Compton scattering of background photons. The mean energy of the secondary photons 
is degraded until it drops below the threshold for pair production, ~ 10 eV for interaction with 
microwave photons. The development of an electro-magnetic cascade, typically over a ~ lOMpc 
distance, converts most of the secondary particles energy to ~ TeV photons. The distance out 
to which > TeV cascade photons may be observed is limited since they may further interact 
with infrared background photons. Current limits on the infrared background radiation energy 
density constrain the mean free path for pair production on infrared photons to lie in the range 
~ 0.6 - 2Gpc for ITeV photons and ~ 20 - 200Mpc for 10 - lOOTeV photons ( [Stecker et al. 1992|) . 

We now estimate the arrival time delays of cascade photons resulting from deflections by the 
IGMF of the protons initiating the cascades and of the pairs produced during the cascades. We 
consider the propagation of protons and secondary particles through an IGMF where a fraction 
r] of intergalactic space is occupied by structures of typical size A and coherent magnetic field B, 
and where the rest of space is filled with a much weaker field. 

First, consider the deflection of a proton with energy E p . After traveling a distance d 
through the IGMF, the proton accumulates a typical deflection angle 8 P ~ (rjd/ X) 1 ' 2 X/ Rl, where 
Rl = Ep/eB is the Larmor radius. This deflection results in a time delay 



where E p = 10 20 E p>2 oeV, d = lOOdiooMpcr? = 0.2rj . 2 , B = 10" n -B_nG and A = 3A 3 Mpc (the 
2/3 factor is due to random field orientations). For IGMF parameters that might be expected 
in models where galactic fields result from the compression of an IGMF produced by turbulence 
induced in the formation of large-scale structure, i] ~ 0.2, A ~ IMpc and B ~ 4 x 10" n G, the time 
delay accumulated by CR protons propagating ~ lOOMpc distance is large enough to reconcile 




(1) 
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the observed CR event rate with the rate of nearby GRBs (r depends sensitively on E p , so that 
the stochastic proton energy loss via pion production results in a broadening of the CR pulse 
over a time comparable to r.) The time delay for cascade photons may be much shorter than 
the CR delay, since the protons lose a significant fraction of their energy over a distance much 
shorter than the ~ lOOMpc distance over which they accumulate their total time delay. Protons 
of initial energy > 2 x 10 20 eV lose 10% of their energy, and therefore initiate electro-magnetic 
cascades carrying ~ 10% of the total cascade radiation energy, over a distance di n it ~ IMpc. 
If the secondary photons and subsequent cascade particles produced in these cascades were to 
propagate rectilinearly to the observer, the time delay of the resulting cascade photons would be 
— di n it/6RL c ~ 0.2i? 2 . 11 day if the burst went off in a high field region, or much shorter if the burst 
happened to go off in a region of low field. If bursts occur inside galaxies, regions of potentially 
much higher field (fiG on a kpc scale), the initial time delay accumulated by 3 x 10 20 eV protons 
as they traverse and leave the host galaxy is only ~ 1 day. 

The pairs produced in a cascade are also subject to deflection by the IGMF, and this can 
lead to much larger photon arrival time delays. The Compton energy loss distance for pairs 
of energy E e = 10 18 , 10 15 eV (corresponding to the high, low energy ends of the cascade), is 
Ajc — 1, O.OlMpc. If these pairs were created in a high field region, their deflection by the 
magnetic field, 6 e ~ eBAic/E e , would delay subsequent cascade photons by ~ 10 3 , 10 5 -B 2 _ 11 yr. 
Thus, the only photons arriving with reasonably short time delays are those produced by branches 
of the cascades where pairs were produced in IGMF voids. (A minimum delay shorter than day 
requires a "void" field < 10~ 15 G). The following characteristics of electro- magnetic cascades 
development ensure, that a significant fraction of the cascade photons are produced by such 
branches, provided that a large fraction of the inter-galactic space is occupied by large scale, 
~ 5Mpc, voids. The distance over which CR protons decelerate, ~ lOMpc, and the mean free 
paths of photons at the high-energy (10 18 — 10 19 eV) part of the cascade, ~ IMpc, are comparable 
to or larger than the typical size of the high field regions. This ensures that a significant fraction 
of the cascades are initiated inside voids, and also implies that, at the high energy end of the 
cascade, a particular cascade branch can "skip" over high field regions without accumulating time 
delay. The short mean free paths of particles at the low energy part of the cascade, Ajc < O.lMpc 
for E e < 10 17 eV, then imply that the development of a branch, whose ~ 10 eV particles are 
created in a void region, is likely to be completed in this region. 

3. Numerical Results 

We now present numerical results for the flux of delayed cascade photons expected in an 
IGMF where a large fraction of inter-galactic space is occupied by large-scale voids, and the 
field strength in the high field regions is high enough to produce a ~ lOOyr delay in the CR 
arrival time. We constructed a Monte Carlo code that propagates high energy nucleons through 
the microwave background and calculates the evolution of the cascades they initiate. All the 
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relevant particle interaction processes are included, and the gyration of charged particles about the 
IGMF is explicitly calculated. Since we are only interested in source distances of order 200Mpc, 
cosmological evolution effects were not included. The cosmic background radiation field was 
chosen to be a superposition of: (i) a blackbody of temperature 2.7°K; (ii) an infrared/optical 
power law distribution at energies e > 0.02eV with density n(e) = 8 x lCP 4 (e/eV)~ 2 cm~ 3 eV _1 , 
corresponding to the lower bound of current estimates of the intensity; (iii) a truncated power law 
distribution at radio energies with density n(e) = 3.5 x 10 _7 (e/eV) -1 ' 75 exp[— e/10 -5 eV]cm~ 3 eV -1 



(e.g., Sironi et al. 1990 ). Our results depend only weakly on the radio background intensity, since 
the cascades are initiated by particles with energies <J 10 19 eV, for which interaction with the 
radio background is only marginally important. The infrared/optical radiation does not influence 
significantly the < lOTeV flux, but may strongly affect the flux at higher energies. This is further 
addressed below. 

We generated several IGMF structures by sampling the sizes of high (low) field regions along 
the line of sight from exponential size distribution with mean 2Mpc(10Mpc). The field in each 
region was randomly oriented with mag nitude B = 4 x 10~ n G in the high field regions and 10 20 G 
in the "void" regions (the "void" field can be as high as ~ 10 -16 G without significantly affecting 
the results shown). The resulting distributions for > 10 20 eV CR arrival time delays are shown 
in Fig. 1 for five different realizations of the IGMF. The source was chosen to be at a distance 
of lOOMpc, and assumed to produce a power law proton spectrum dN/dE oc E~ 2 in the range 
10 20 - 10 21 eV. The corresponding photon fluences above ITeV are shown in Fig. 2. The fluence is 
normalized for a burst producing 10 51 ergs in > 10 20 eV protons, as would be required to account 
for the observed CR flux. The cascade radiation spectrum is determined by the development 
of the low energy end of the cascade, and is therefore independent of the proton production 
spectrum. For the short time delays of interest here, the low energy development of the relevant 
cascade branches always occurs in void regions. Thus, the spectrum is similar to that obtained 
in the absence of an IGMF and is time independent. The integral photon number flux (above 
energy E) is approximately given by /(> E) oc E~ - 8 for 3GeV < E < 200GeV, /(> E) oc E^ 1 - 
for 200GeV < E < lOTeV, and /(> E) oc E~ 1A for lOTeV < E < 70TeV. 

The spread in the typical CR time delay obtained for different field structure realizations 
is small, since this delay is accumulated over a lOOMpc propagation. The spread in the photon 
fluence obtained on a time scale of days is larger, since it is determined by the IGMF structure near 
the source, at distances < 20Mpc, where the development of most of the relevant electro-magnetic 
cascades take place: the highest fluence was obtained for a case where the first 20Mpc were free of 
high field regions, and the lowest for a case where 50% of the first 20Mpc were occupied by two 
~ 5Mpc high field regions. This reflects the results of a large set of Monte-Carlo calculations, 
which show that for B ~ 4 x 10~ n G and void regions size in the range 5-10Mpc, the fluence 
increases from ~ 10~ 7 cm~ 2 to ~ 10~ 6 cm~ 2 as the fraction of the volume occupied by high field 
regions within 20Mpc from the source decreases from 0.4 to 0.2. The presence of a host galaxy 
near the burst is expected to introduce a minimum delay of order lday in the cascade photon 
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arrival times, and is not expected to significantly affect the fluence on longer time scales. This is 
demonstrated in Fig. 2 by the two Monte-Carlo calculations, shown with the heavy lines, one of 
which includes a 1/iG magnetic field over the first lkpc propagation distance to simulate the effect 
of a host galaxy magnetic field. It should be noted that the host galaxy delay can be considerably 
shorter if the burst produces a significant amount of energy in protons with energies ^ 10 21 eV, for 
which a host galaxy would introduce only a ~ lh delay. 

Photons arriving with ^ lday delay are produced by cascades initiated by pions created over 
the first lOMpc propagation distance of the protons. Most of these pions are produced by protons 
with energy > 2 x 10 20 eV, which lose ~ 1/2 their energy over this distance. The typical > ITeV 
fluence obtained on a day time scale, ~ 10~ 6 cm~ 2 (~ 10~ 5 erg cm -2 ), corresponds to ~ 2% of the 
energy produced by the burst as > 2 x 10 20 eV protons. 



4. Conclusions 

Although a scenario in which CRs above 10 20 eV are produced by cosmological GRBs requires 
the arrival time of CRs to be delayed with respect to the 7-rays by more than ~ 50yr, we have 
shown that the delay of secondary, 0.01 — lOOTeV cascade photons may be much smaller. A 
short delay is possible provided that a large fraction of the inter-galactic medium is occupied by 
large-scale magnetic "voids", regions of size ^ 5Mpc and field weaker than 10~ 15 G. For a field 
strength of ~ 4 x 10~ n G in the high-field inter-galactic regions, which would account for the 
observed galactic fields if it were frozen in the protogalactic plasma, the delay of CRs produced 
by a burst at a distance of lOOMpc is ~ lOOyr. At the same time, the fluence of secondary 
> ITeV photons on a lday time scale would be ~ 10 _6 cm~ 2 , provided that ~ 80% of the 20Mpc 
region around the source is occupied by magnetic "voids" (The fluence is inversely proportional 
to the burst distance squared, since the photon time delay is independent of the burst distance). 
The integral photon number flux in the energy range lOGeV < E < lOTeV is approximately 
/(> E) (x E- 1 - . The flux at higher energies is very sensitive to the infrared background intensity. 
If the intensity is near the lower bound of current estimates, the flux extends beyond lOTeV, 
approximately as /(> E) cx E~ 1A . If the intensity is close to its current upper bound, it would 
completely suppress the > lOTeV flux from distances > lOOMpc. 

A 3<7 detection of > ITeV photons by current high-energy 7-ray experiments requires a fluence 
~ 10~ 6 cm -2 y/tdayEmln TeV , where td ay is the observation time measured in days and -EVmn.TeV is 



the detector threshold energy in TeV (see Alexandreas et al. 1991 for CYGNUS, Funk et al. 1995 



for HEGRA, |McKay et al. 1993; for CASA-MIA, Cawley fc Weeks 1995| for Whipple). This fluence 



is close to that expected from a burst at a distance of ~ lOOMpc. However, in a cosmological 
model, the rate at which GRBs occur in a lOOMpc sphere around us is low, ~ O.lyr -1 . A factor 
of 10 increase in the sensitivity of TeV detectors, as expected for example in the near future in 
the Whipple and HEGRA observatories, would allow the detection of the delayed flux from bursts 
occuring at distances up to ~ 300Mpc. The rate of such bursts is ~ 2yr _1 . 
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Lower threshold energy, space-based detectors such as EGRET may also detect the delayed 
flux. At lOGeV, EGRET has an effective area of ~ 10 3 cm 2 . Thus, for the > lOGeV fluence 
expected from a burst at a lOOMpc distance, ~ 10~ 4 cm~ 2 in one day, the probability that EGRET 
detects a > lOGeV photon is ~ 0.1. This probability, although not negligible, is small. Therefore, 
detection of the delayed flux would probably require next-generation GeV instruments, such as 



GL AST ( |Bloom 1995|) , that are expected to have order of magnitude better sensitivity [It should 



be noted, however, that EGRET has detected a 18GeV photon from the direction of one of the 



strongest BATSE bursts (second in fluence), with ~ 1.5hours delay ( Hurley et al. 1994 )] 
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Fig. 1. — The distribution of arrival time delays of > 10 eV protons produced by a burst at a 
distance of lOOMpc, for various IGMF structures where a fraction ~ 0.2 of the inter-galactic space 
is occupied by ~ IMpc regions of coherent 4 x 10~ n G field. 
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Fig. 2. — The fluence of > ITeV secondary photons, following a GRB producing 10 51 erg as 
> 10 20 eV protons at a distance of lOOMpc. The IGMF structures are the same as those used in 
Fig. 1, and similar line-styles correspond to the same field structure in both figures. The heavy 
dotted line gives the fluence obtained when a 1/j.G host-galaxy magnetic field is added to the field 
structure corresponding to the heavy solid line. 



